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Abstract 



t4-H ' For sedimenting colloidal hard spheres, the propagation and broadening of the crystal-fluid in- 

O ' 

. ' terface is studied by Brownian dynamics computer simulations of an initially homogeneous sample. 



Two different types of interface broadenings are observed: the first occurs during growth and is 



i-rt '. correlated with the interface velocity, the second is concomitant with the splitting of the crystal- 

Q ' fluid interface into the crystal-amorphous and amorphous- liquid interfaces. The latter width is 



strongly peaked as a function of the gravitational driving strength with a huge amplitude relative 
to its equilibrium counterpart. 
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I. INTRODUCTION 

Sedimentation of colloidal particles in a liquid is a widespread phenomenon which governs 

n 

the formation of river sediments close to salty sea water [1], controls the function of red blood 
cells a. a„d ,s tech„olog,cally used to separate d,ffere„t sorts of part.cles by ceutrifugation 
p[. Gravity is also used to compactify colloidal samples. In a finite container, the actual 
settling process reaches a sediment at ion- diffusion equilibrium, which is characterized by a 
static colloidal density profile. Much of our particle-resolved knowledge of dense suspensions 
stems from model hard sphere colloidal dispersions made up of sterically-stabilized particles. 
In equilibrium, it has been shown that density profiles contain the isothermal equation of 
state of the hard sphere system [4, l5|. At high gravitational strengths (or Peclet numbers) 
the density of hard spheres at the bottom of suspension increases above the lower limit of the 
solid-liquid coexistence (f) = 0.492 p, and a crystallization in the bottom layers occurs [7H9|. 
The number of crystalline layers is controlled by gravity and by the sedimentation velocity 
of particles, and the overall particle density per area. While the equilibrium properties of 
crystal sediments are well-explored by now, the dynamics and relaxation towards equilibrium 
is by far less understood [10|, lll|. In particular, it is known that crystallization occurs 



under gravity at different conditions than in the bulk 12|, [iSl]. The technique of colloidal 
templating has been used to steer colloidal crystallization layer-by- layer in gravity [ij, [l5| 



and the kinetics of crystalline defects [16 



llTj has been debated for sediments. 



Using real-space confocal microscopy techniques, DuUens and coworkers [18| have recently 
measured a dynamical broadening of the hard sphere solid-fluid interface during sedimenta- 
tion. This points to the essential role of crystal nucleation in the supersaturated fluid and 
the subsequent built-in of "crystal packages" into the interface. In a very slowly growing 
interface, and at low Peclet numbers, a broadening of the fluid-solid interfacial width relative 
to its equilibrium value was found. In this letter, we address the broadening of the solid-fluid 
hard sphere interface in the sedimentation process by extensive Brownian dynamics com- 
puter simulations. We find two different types of broadenings of the fluid-crystal interface: 
the first occurs during growth and is correlated with the interface velocity. This type of 
broadening is not driven by gravity as it also shows up in the gravity- free case [l9[. The 
second type of interface broadening occurs before the splitting of the crystal- amorphous 
interface from the amorphous-liquid interface. The separated interface includes a region 




FIG. 1. Laterally averaged packing fractions 4'{z,t) (thin pink lines) versus reduced height z/a 
for three different times (a-c): t = 20TB,40TB,100rB. The parameters are Pe = 0.5, r]A = 45.4, 
and 0=0.3. Also shown are the interfacial profiles nc{z,t) for the "crystalline" particles (thick red 
line) and ni{z,t) for the " liquid- like" particles (dashed blue line). The interface position zo{t) is 
determined from the crossing point of the thick red and dashed blue lines. 



which is structurally disordered and dynamically arrested. The latter width behaves non- 
monotonic in the gravitational driving strength (or Peclet number Pe) and exhibits a huge 
peak relative to its equilibrium counterpart. While the first type of broadening has not yet 
accessed in confocal microscopy, the latter is consistent with the measurements of DuIIens 
et al [18| . However, their measurements were only performed for low Peclet numbers before 
the splitting of two interfaces happen. 
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FIG. 2. Reduced interface position zo(t)/a (a), reduced interfacial width w{t)/a (b), and reduced 
propagation velocity v{t)TB/(y (c) as a function of reduced time t/rB- Full blue lines with circles 
are for r]A = 45.4 and (f) = 0.1, full red lines are for tja = 45.4 and (f) = 0.3, dashed pink lines are 
for r]A = 90.8 and (j) = 0.3. Full black lines with squares are simulation data for the gravity- free 
{Pe = 0) suspension with a volume fraction </> = 0.51 and r]A = 46.27. 

II. SIMULATION AND ANALYSIS 



In our Brownian dynamics computer simulations, we use an adapted code for hard spheres 

20|] of diameter a, where the short-time infinite-dilution diffusion coefficient Dq sets the 

Brownian time scale t-q = a"^ / Dq. A time step of At = O.OOItb was used in integrating 

the stochastic equations of motion. The simulation box contains A^ hard spheres and has a 

rectangular shape with dimensions L^ = 40.8cr, Ly = 43. 2o", and L^ varied between 54a" and 



^3 



24:0a. Various systems with a packing fraction = ^^ ^ ^ j^ in the range 0.1 < < 0.45 
were investigated. A further important system parameter is the gravitational load, or the 
surface (or areal) density of particles pa = N/L^Ly, which we scale to 

6L 

7]A = Na^L^Ly = 0^ (1) 

na 

Thus different combinations of and L^ result in the same t]a- Simulations were carried for 
r]A = 45.4 with A^ = 80 000 particles and for r]A = 90.8 with A^ = 160 000 particles. The 
quantity rjA is chosen large enough to produce a crystal-fluid interface at a given gravitational 
strength Pe. The gravity acceleration g points along the — z-direction and has a relative 
strength 

Pe = mga/ {2kBT) (2) 

where ksT denotes the thermal energy and m the buoyant mass of the colloidal particles 



2l| . Periodic boundary conditions are employed in the x and y directions, while two hard 



walls are placed at 2; = and z = L^. Next to the wall at 2; = we place a triangular 
layer of fixed spheres with a lattice constant a = 1.133cr which acts as an initial template 
for crystal growth 22|. Without a template the crystallization happens in a few bottom 
layers containing small grains, large defects and fault stackings. This strongly suppresses 
the formation and following upward propagation of a single-phase crystalline front along the 
sediment. 

All simulations were started from an initial configuration with a homogeneous distribu- 
tion of colloids in the simulation box except the template particles in the seed layer next 
to the bottom wall. This mimics an initially stirred solution to which gravity is applied 
instantaneously [11]. A total simulation time of 500rB was accessed during our simulations. 
Within this time each particle can sediment about lOOOaPe, a distance at least 5 times 
larger than the characteristic sedimentation length h ^ H x 0.55/0. Here H is the hight 
of the sedimented segment above which the packing fraction of suspension is < 0.5. We 
averaged our final stage and time-independent results over 20 different initial configurations 
in order to improve the statistics. 

For subsequent times t, we calculated the laterally averaged one-particle packing fraction 
which is defined as 

^i^^''^) = qJ~J~ J J dx dy p{x,y,z,t) (3) 

5 










Pe =0.5 


(a) 


1 1 1 J 1 1 1 1 1 i i lull 


lim 


J' " 

11 


ill 


w^M«MVu 


Pe=1.6 


(b) 


1 

M 


ll^^ 


X 


Pe =3.0 


(c) 




\, 







10 



20 
z/a 



30 



40 



FIG. 3. Same as Figure [H but now for a fixed large time t = ti = 500tb and three different Peclet 
numbers: (a) Pe = 0.5, (b) Pe = 1.6, (c) Pe = 3.0 at r]A = 45.4 and (j) = 0.1. 



and is resolved along the vertical 2;-coordinate. Here p{x, y, z, t) is the local one-body density 
of particles at a given time t. We further identified "crystalline" particles with a crystal- 
like surrounding according to a commonly used criterion 23|: the local orient at ional-order 
parameter qQ{i) is calculated for each particle i 2J, |25|]. When two particles i and j are 
separated by r < 1.3a, we associate a crystalline bond to these particles, if (jQ^i) ■%{]) > 0.5. 
A particle which has at least 8 of these bonds is considered to be crystalline. All other 
particles are identified as a " liquid-like" particles. The corresponding local packing fractions 
of crystalline particles nc{z,t) and liquid particles ni{z,t) are also calculated and give rise 
to a Qq interface, an example of which is given in Figure [H We now define the interface 
position Zo{t) by the implicit condition nc{zo{t),t) = ni{zo{t),t), i.e. by the position where 



the fraction of crystalline-particles equals that of the liquid-like particles. The interfacial 
width w{t) is defined as the inverse of the maximal slope of the tIc profiles 



'dnc{z,t)\ 
dz 



-W=(^^^^§e^l UoW (4) 



We finally define the propagation velocity of the solid-fluid interface as v{t) = dzo{t)/dt. 

III. MECHANISM 1: KINETIC BROADENING 

Results for the interfacial profiles at three different times are shown in Figure [U Clearly 
the crystal-fluid interface is propagating from the bottom of the container (left side in the 
figure) into the fluid. The lowering of peak heights in the density profile of Figure [Th in the 
region between z = lAa and z = 27a originates from the stacking faults |9|, |26| and defects 
introduced by a small amount of rhcp and hep crystalline grains. These defects tend to 
dissipate into fee crystalline layers during longer tirne simulations. The system size also plays 
a crucial role in the dissipation of stacking defects [27]. When the crystal grows, the interface 
is getting broader. This is clearly illustrated in Figure |2] where the interface position ZQ{t), 
its width w{t) and its propagation velocity v{t) are shown simultaneously. The interface 
velocity is non-monotonic in time: it first gets accelerated, reaches a maximum and decays 
to zero. The initial acceleration has to do with an "induction time" to get particles into 
positions which are structurally favorable for subsequent crystal growth. The final slowing- 
down is due to the approach to equilibrium where ideally a finite height of the interface is 
reached. Concomitantly with the non-monotonic behavior of the interface velocity, there is 
a broadening of the interface: as a function of time, it first broadens and then shrinks again. 
The saturation value at large times, however, is significantly larger than the initial width. A 
direct comparison of Figure Eb andEt reveals that the maximum of the propagation velocity 
coincides with the maximal interfacial width demonstrating that these two phenomena are 
correlated. The effect is stable for different overall densities rj^. An increased rj^ just retards 
the occurrence of the velocity and width maximum such that they occur at larger simulation 
times. A similar retardation effect takes place in systems with low packing fraction at a 
fixed areal density t]a, compare solid lines with and without symbols in Figure [2l 

As a first type of interfacial widening, we therefore identify a kinetic broadening of the 
interface: upon growth, the qQ-mteiiace becomes less steep, see thick lines in Figure [TJd. 



This has to do with the fact that the system has less time to structurally arrange when the 
interface propagates quickly. An additional simulation for zero gravity revealed that the 
kinetic broadening effect exists also at zero Peclet number [28|. However, in the gravity- 
free case the front growth velocity has no maximum. As a reference, simulation data for 
the front position, interfacial width and the propagation velocity for a coexisting fluid and 



23 



29 



30| in a gravity-free system are included into Figure [2] as full lines with 



crystal 
squares. 



IV. MECHANISM 2: INTERFACE SPLITTING 

We now address the crystal-fluid interface in the final state of the sediment after a long 
time ti = SOOtb and discuss its properties as a function of the Peclet number (at fixed 
rjA)- In Figure 131 the laterally-averaged packing fraction profiles are shown at t/ for three 
different Peclet numbers Pe = 0.5, 1.6, 3.0. Interestingly, the position of the interface (as 
visualized by the crossing of the red and blue lines) behaves non-monotonic in the Peclet 
number, see also Figure HJ which contradicts the expectation based on equilibrium argu- 
ments. This expectation is based on a discussion of the two limiting cases Pe — )■ and 
Pe — 7- cxD in equilibrium: for zero gravity, the density profile is homogeneously distributed 
over the full simulation box (apart from local density correlations at the hard wall) with a 
low bulk volume fraction. Hence the solid-fluid interface (if at all) is close to the bottom 
wall. In the opposite limit Pe — )■ oo, there is a finite number of layers with a closed packed 
density which then abruptly drops to a vacuum at higher Pe at a height of about 30cr. In- 
terpolating between these two limiting cases, in equilibrium, the interface position increases 
with Peclet number Pe. Starting from Pe ~ 0.8 the single crystalline-liquid interface splits 
into two separate interfaces, the crystalline- amorphous and the amorphous-liquid interfaces. 
When the gravitational load increases, the amount of amorphous glass grows partly from the 
sedimented liquid particles, and partly from the collapsed crystalline layers. Higher loads 
completely destroy the crystal leaving only a few intact layers. Besides of the gg-interface 
position, a position of the amorphous-liquid interface can be defined at the height where 



the coarse-grained density (averaged over the oscillations) 3l|] equals the bulk coexisting 
fiuid volume fraction of 0.492. As opposed to the nonmonotonic variation of the gg-interface 
position, the position of the amorphous-liquid interface is monotonic, see Figure HI While 
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FIG. 4. (a) Position ZQ{ti)/a of the gg-interface (full lines), and of the amorphous-liquid interface 
(dashed lines) versus Peclet number Pe for two areal densities tja = 45.4 (a) and r]A = 90.8 (b) at 
a fixed large time t = ti = 500tb- The interface splitting is clearly visible. Black squares at Pe = 2 
correspond to the amorphous-liquid and (/g-interfaces analyzed in Figure [5l 



the positions of the qQ-mteriace and the amorphous-hquid interface coincide for small Peclet 










FIG. 5. Simulation results for the sediment at t = t/ = 500tb for system parameters Pe=2, 
4> = 0.3 and rjA = 90.8. (a) Schematic picture showing the crystalline, amorphous and fluid zones 
in the sediment, (b) Laterally averaged packing fraction 0(z) (thin pink lines), interfacial profile 
nc{z) for the "crystalline" particles (thick red line), interfacial profile n;(z) for the " liquid- like" 
particles (dashed blue line) versus reduced height z/a. (c) Mean square displacement A(z) of the 
particles along the sediment obtained within a time window of At = lOr^. The positions of the 
crystal-glass and glass-fluid interfaces correspond to the black squares in Figure IDa. 



number, there is a splitting above a threshold leading to the non-monotonic behavior. The 
non-monotonicity can be explained by the formation of a dynamically arrested region be- 
tween the splitted q^ and amorphous-liquid interfaces schematically shown in Figure [5^ for 
a system with cf) = 0.3, rjA = 90.8 and Pe = 2. Within the glassy amorphous region the 
particles are structurally disordered, see Figure Eb, and dynamically caged such that they 
cannot find suitable surrounding to nucleate into a big, and layered, crystal. This is doc- 
umented by the smallness of the mean square displacement of particles in the amorphous 
layer, see Figure ^. The amorphous dynamically arrested part stops the propagation of 
the gg-interface. Strong gravity acts therefore like a fast and deep compression, a situation 

10 



which favours glass formation in general. The disordered region does not seem to be com- 
posed of polycrystalline material as this would have resulted in a finite fraction of crystalline 
particles. As shown in Figure S^ andSb, the Peclet number threshold for splitting decreases 
for increasing packing fraction (f) at fixed r^^, but does not depend on the imposed overall 
density tja for fixed 0. Closely above the Peclet number threshold at which interface split- 
ting occurs, the width w exhibits a strong peak, see Figure [61 In fact, the width increases 
by an order of magnitude relative to its low-Peclet-number value and then decreases again. 
This surprising huge broadening points to an extremely broad interface between the initial 
crystalline layers and the subsequently formed disordered material. The subsequent decrease 
for higher Peclet number is qualitatively understood by considering the limit Pe — )■ oo: if 
everything is dominated by strong gravity, the interface then is expected to sharpen due to 
enforced crystallization. 



We now compare our results to the real-space experiment by Dullens and coworkers 18|. 
The broadening observed in Ref. [l8| was measured for long times at relatively small Peclet 
numbers (smaller than 0.8) and for larger gravitational loads tja than considered in our work. 
The interfacial velocity was small but nonzero. This implies that the broadening observed 
in [l8| is consistent with our findings at long times for low Pe, before the interface splitting 
occurs. It would be interesting to do further experiments at higher Pe to verify the interface 
splitting and the huge interface broadening predicted by our simulations. 

V. CONCLUSION 

In conclusion, we have studied the width of a crystal-fluid interface in a sedimenting 
suspension of hard spheres by Brownian dynamics computer simulations. Two qualitatively 
different types of broadening are observed. The first type of interfacial broadening is purely 
kinetic. It is not triggered by gravity but just correlates with the interface velocity. The 
second type of broadening is almost static and is huge if the structural crystallinity interface 
and the amorphous-liquid interface split. Between these two interfaces, an amorphous dense 
sediment is formed which exhibits a wide structural interface with the lower crystalline par t 



of the sample. This prediction can in princip 
It might also be important for granulates 



3e verified by real-space experiments 18 



23|. 



crystallization in complex plasma 33 1 



32( 1 where the Peclet number is high, and for 
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FIG. 6. Interfacial width w{ti)/a versus Peclet number Pe at a fixed large time t = ti = 500tb 
for two different packing fractions (/> = 0.1 (line with symbols) and (p = 0.3 (full line) at fixed areal 
density rjA = 45.4, and for packing fraction = 0.3 and rjA = 90.8 (dashed line). 



Real experimental samples possess an intrinsic polydispersity. Furthermore there are 
solvent-mediated hydrodynamic interactions. Most of the influence of hydrodynamic inter- 
actions can be captures by an effective short-time diffusion constant Dn. Nevertheless the 



influence of polydispersity 23|] and hydrodynamic interactions 11 



34 



35| on crystal growth 



remains for future exploration. Furthermore new effects are expected for attractive systems 
36( 1 as witnessed by different dynamically arrested states (attractive glasses) 37| which could 
impede crystallization during sedimentation. The latter case is realized for colloid-polymer 
mixtures. 

It is worth to mention that some of the present flndings on the interface broadenings 
may be characteristics of fee {111} growth from a patterned template with triangular lattice 



12 



considered in current study. It will be useful to check whether a similar broadenings takes 



place in the {100} growth, a case considered in Refs. [29l. l30j| for gravitation- free hard sphere 
systems. 

In general, our findings can be used to steer the thickness of crystalline and amorphous 
layers as well as their interfacial structure by sedimentation. This may help to open the way 
for designing smart colloidal materials with new optical and rheological properties. 
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